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Abstract — We transformed radial velocities compiled from more than 1400 published sources, including the Geneva- 
Copenhagen survey of the solar neighborhood (CORAVEL-CfA), into a uniform system based on the radial velocities of 
854 standard stars in our list. This enabled us to calculate the average weighted radial velocities for more than 25 000 
HIPPARCOS stars located in the local Galactic spiral arm (Orion arm) with a median error of ±1 km s -1 . We use 
these radial velocities together with the starscoordinates, parallaxes, and proper motions to determine their Galactic 
coordinates and space velocities. These quantities, along with other parameters of the stars, are available from the 
continuously updated Orion Spiral Arm CAtalogue (OSACA) and the associated database. We perform a kinematic 
analysis of the stars by applying an Ogorodnikov-Milne model to the OSACA data. The kinematics of the nearest single 
and multiple main-sequence stars differ substantially. We used distant (r ~ 0.2 kpc) stars of mixed spectral composition 
to estimate the angular velocity of the Galactic rotation, lu = —25.7 ± 1.2 km s -1 kpc -1 , and the vertex deviation, 
I = 13° ± 2°, and detect a negative K effect. This negative K effect is most conspicuous in the motion of A0-A5 giants, 
and is equal to K = —13.1 ± 2.0 km s -1 kpc -1 . 
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1. INTRODUCTION 

Radial velocities of several tens of thousands of 
stars have been measured in recent years with preci- 
sions of better than ± 1 km s , making joint analysis 
of stellar proper motions and radial velocities of con- 
siderable importance for studies of stellar kinematics. 
An appreciable fraction of these high-precision radial 
velocities were measured using the CORAVEL pho- 
toelectric cross-correlation spectrometers and digital 
spectrometers of the Harvard-Smithsonian Center for 
Astrophysics (CfA) and can be found in the catalog of 
Nordstrom et al. [1] (hereafter, the CORAVEL-CfA 
catalog). These authors and Gontcharov [2] showed 
that the combined use of proper motions and radial 
velocities, interstellar extinction estimates, multicolor 
photometry, and binary parameters makes it possible 
to draw conclusions about the absolute magnitudes, 
ages, masses, and Galactic orbits of the stars and the 
evolution of Galactic structures over millions of years, 
and to perform other stellar-kinematic studies. In or- 
der to systematically collect, analyze, and reduce the 
corresponding observational material, we have created 



and have been continuously updating our database for 
stars with known coordinates, parallaxes, proper mo- 
tions, and radial velocities. 

Although the mutual alignment and density of spi- 
ral arms in the solar neighborhood require a special 
study, intrinsically bright HIPPARCOS stars, which 
are tracers of the Galactic spiral arms, enable us to 
estimate fairly confidently the position and size of the 
local spiral arm. The HIPPARCOS catalog is complete 
to about 7.3 m in the band of the Milky Way, i.e., for 
stars with absolute magnitudes M < —2 m located out 
to 600 pc, given that, as a rule, the extinction is lower 
than 0.5 m in this region. Moreover, stars with absolute 
magnitudes M < — 4 m belonging to neighboring spiral 
arms are represented almost completely in this catalog. 
The distribution of intrinsically bright stars over the 
Galactic octants and in heliocentric distance present- 
ed in Table 1 indicates that, first, nearby intrinsically 
bright stars are distributed almost uniformly over the 
Galactic octants; second, such stars are not distribut- 
ed uniformly in Galactic longitude at greater distances; 
and, third, the distribution changes abruptly at a helio- 
centric distance of more than 700 pc. Thus within 400 



pc from the Sun, only the Scorpio. Sagittarius associa- 
tion can be identified near a Galactic longitude of about 
340.. At distances between 400 to 700 pc, intrinsical- 
ly bright stars are mostly concentrated near latitudes 
of 90° and 225°, possibly indicating the directions to 
which the local spiral arm extends, as well as the width 
of the arm, which is about ±500 pc in the directions of 
the Galactic center and anticenter. The predominant 
longitudes change beyond 800 pc: the regions of en- 
hanced density of intrinsically bright stars near 110° 
and 290° correspond to the well-known Perseus and 
Sagittarius arms, respectively. We can thus refer to 
the region within 500 pc of the Sun as the portion of 
the local Orion arm that surrounds us. 

Currently, radial velocities are known and paral- 
laxes are relatively accurate primarily for stars within 
500 pc of the Sun, i.e., within the local Galactic spiral 
arm (Orion arm). For this reason, we called the catalog 
based on our database the Orion Spiral Arm CAta- 
log.OSACA Jhttp:// www.geocities.com /orionspiral / ) . 
The database, with its numerous tables, relations, and 
queries is a sort of "factory", whereas the final product 
that can be distributed is the OSACA catalog, which 
consists of several tables that are periodically compiled 
from the database. 

The goal of this paper is to describe the catalog and 
to analyze the most interesting examples of its use for 
studies of the kinematics of stars in the solar neighbor- 
hood. One such example is analyzing the effect of an 
adopted model for the binarity of the nearest stars on 
their kinematic parameters. As another example, we 
considered relatively distant A0.A5 stars exhibiting a 
negative K effect. We performed our kinematic anal- 
ysis using a Ogorodnikov-Milne model. We determined 
the Schwarzschild ellipsoid to analyze the residual ve- 
locities of the stars and construct two-dimensional dis- 
tributions in the U, V plane. 



2. CONTENT AND FORMAT 
OF THE OSACA CATALOG 

The core of the OSACA catalog are the six quanti- 
ties that describe the position and velocity of a star in 
space: a, 5, tt, fi a , fi>s, and V r . These quantities are used 
to calculate the coordinates X, Y, Z and space veloci- 
ties U, V, W [3] in the standard Galactic Cartesian co- 
ordinate system, where X increases toward the Galactic 
center, Y in the direction of Galactic rotation, and Z in 
the direction of the North Galactic Pole. The OSACA 
catalog is continuously updated with new published 
data. Every substantial change in the input data is fol- 
lowed by a recalculation of the dependent quantities; 
for example, the velocity components U, V, and W are 
recomputed after a refinement of the radial velocities 
in the database. 

Since positions and velocities are known sufficient- 
ly precisely only for stars of the HIPPARCOS catalog 



[4], the OSACA catalog currently includes only these 
stars, although we collect data for other stars in our 
database as well. The proper motions of the HIPPAR- 
COS stars considered have median errors in proper mo- 
tion of ±0.0007" yr _1 , in right ascension and declina- 
tion of ±0.0006" yr _1 , and in the tangential velocity of 
±0.0009" yr" 1 . The proper motions of some of the stars 
considered have low accuracy due to their binarity. To 
refine these proper motions, we used a combination of 
coordinate and proper-motion data from the HIPPAR- 
COS catalog and from ground-based catalogs, in accor- 
dance with the method of Gontcharov and Kiyaeva [5]. 
The parallaxes, equatorial coordinates, and proper mo- 
tions of the stars considered are updated comparatively 
infrequently, when new data whose accuracy are supe- 
rior to the HIPPARCOS data become available. Thus, 
the main data that are collected and refined are the 
stellar radial velocities. 

The catalog currently includes the main table, a ta- 
ble of references to published radial velocities (one line 
per star) , a table of references to published radial ve- 
locities (one line per publication), and a table of cross 
identifications of the stars. Table 2 gives the format 
of the main table of the current release of the catalog. 
Figure 1 shows the distribution of the stars in absolute 
magnitude and heliocentric distance for single (solid) 
and multiple (dashed) stars. The distribution of the 
stars in the absolute magnitude. B.V color index dia- 
gram (Fig. 2) resembles the corresponding distribution 
for all stars of the HIPPARCOS catalog. all spectral 
types and luminosity classes are represented. 



3. CALCULATION OF THE RADIAL 
VELOCITIES 

The WEB catalog [8] contains about 14 300 HIP- 
PARCOS stars with precise parallaxes and radial ve- 
locities known to better than ±5 km s , The cata- 
log of Barbier-Brossat and Figon [9] contains about 16 
000 such stars. Our catalog contains more than 28 000 
such stars, reflecting the large number of high-precision 
radial- velocity measurements obtained in recent years, 
including the measurements for more than 10 000 stars 
reported in the CORAVEL-CfA catalog, as well as our 
reduction of all the radial-velocity catalogs employed 
to a single system. 

We calculated the weighted mean radial velocities 
using data from more than 1400 publications. The full 
list of these sources can be found on the OSACA web 
page (http://www.geocities.com/orionspiral/). Exam- 
ples of extensive catalogs that we employed include 
[10.21] and all the sources used to produce the cata- 
logs [8, 9, 22]. All the sources of the radial velocities 
can be subdivided into two parts: (1) about one hun- 
dred major source catalogs, each of which usually has 
its own system of radial velocities, and (2) over one 
thousand publications dedicated to individual stars or 



small groups of stars. In the former case, the system- 
atic errors for a particular catalog can be revealed by 
comparing it to some standard or other catalog, where- 
as such a comparison is impossible in the latter case 
due to the small number of stars. If individual obser- 
vations of stars were made with the same instrument 
as the measurements reported in some major catalog, 
we can hope that the radial-velocity system will be the 
same for these observations, so that the systematic cor- 
rections assigned to the instrument in question can be 
applied to that individual star. 

The problem of radial- velocity standards has a long 
history, and has become especially topical in recent 
years (see, e.g., [23.27]). At present, only 181 stars ap- 
proved by the International Astronomical Union (IAU) 
as radial-velocity standards have been confirmed to 
have constant radial velocities. The HD and HIPPAR- 
COS numbers, radial velocities, V magnitudes and ab- 
solute magnitudes, and B.V color indices of these stars 
are given in the table posted on the OSACA web page; 
in these cases, the "source "column indicates "IAU". 

However, the need to use certain stars as radial- 
velocity standards has been questioned in recent years 
due to the emergence of astrometric methods of mea- 
suring radial velocities [16] and the fundamental pos- 
sibility of calibrating an instrument by observing the 
radial velocities of Solar-system bodies. 

Nidever et al. [27] (which we will refer to as the 
Keck catalog) shows that: (1) the use of standard stars 
and Solar-system bodies yields virtually the same ra- 
dial velocities to within ±0.5 km s _1 , (2) it is possi- 
ble to create a list of several hundred standard stars 
for determining radial velocities with errors no greater 
than ±0.5 km s -1 , which would reasonably represent 
the entire variety of spectral types, luminosity classes, 
and other stellar parameters, and (3) determining the 
absolute radial velocities to better than ±0.5 km s" 1 
requires allowance for numerous effects that thus far 
lack a sufficient theoretical basis (to say nothing of the 
procedures for taking them into account in practice!). 

Stellar kinematic studies require not so much ra- 
dial velocities with very high precision as a homoge- 
neous data set of radial velocities for a large number 
of various stars. A precision of 1 km s _1 can be con- 
sidered to be acceptable for such studies. A velocity of 
1 km s _1 corresponds to 1 pc Myr -1 . This precision 
enables the calculation of Galactic orbits and the anal- 
ysis of the evolution of Galactic structures and groups 
of stars within 500 pc of the Sun over several tens of 
million of years in the past and future. Moreover, with- 
in a large volume of space, a precision of 1 km s" 1 is 
comparable to the median accuracy of stellar tangen- 
tial velocities calculated from proper motions, namely, 
0.001 " yr" 1 ; within 200 pc of the Sun, the tangential 
velocities are, on average, known better than the cor- 
responding radial velocities (for about two-thirds of all 
the stars considered), whereas the situation is reversed 
for more distant stars (about one-third of all the stars) . 



The Keck catalog includes a list of 782 stars pro- 
posed for use as standard stars for determinations of 
radial velocities. These stars have stable relative radi- 
al velocities, and represent different types of stars; the 
radial velocities of IAU standard stars in the Keck cat- 
alog agree with the standard radial velocities to within 
no worse than ±0.5 km s _1 . Although some of these 
stars exhibited radial-velocity variations after the re- 
lease of the Keck catalog, the overwhelming majori- 
ty can be adopted as supplementary standards for the 
determination of radial velocities with a precision of 
±(0.5-1) km s _1 . The parameters of 673 such stars 
with "Keck"indicated in the "source "column are given 
in the table posted on the OSACA web page. 

The full list of 854 standard stars that we used 
to transform more than 50 000 published radial ve- 
locities to a single standard system that we refer to 
as the "IAU-Keck" system does not include some stars 
in the initial IAU and Keck lists that have exhibited 
appreciable radial- velocity variations in recent years. 
Figure 3 shows the distribution of these stars on the 
(B — V) — My color index. absolute magnitude diagram. 
It is evident that types of stars that are not represented 
or under- represented include white dwarfs, subdwarfs, 
subgiants, and A5 to F5 main-sequence stars. 

The transformation of a published radial-velocity 
catalog into the standard system involves the calcu- 
lation of "catalog minus standard "differences for stars 
in common followed by the identification and approx- 
imation of systematic dependences of these differences 
on various stellar parameters allowing the systematic 
differences identified to be eliminated for all stars. 

As an example, let us consider the systematic er- 
rors that we revealed for the CORAVEL-CfA catalog. 
Figure 4 shows the dependences of the radial-velocity 
differences in the sense "CORAVEL-CfA minus stan- 
dard "for 484 stars in common as a function of radial 
velocity, B.V color index, and absolute magnitude (we 
also found the dependences on effective temperature, 
Hp index, and metallicity Fe/H). The left- and right- 
hand plots show the data before and after elimination 
of the following dependences: 

AV r = -0.8398 * (B —V) + 0.2531, 

AV r = 0.000008 * V r 2 + 0.0012 * V r - 0.0142, 

where SV r is the difference of the radial velocities for 
stars in common in the sense CORAVEL-CfA minus 
standard, B.V is the color index, and V r is the radi- 
al velocity. It is obvious from these plots that elimi- 
nating the dependences on radial velocity and B — V 
color index also eliminated the dependence on abso- 
lute magnitude (and on other parameters as well); no 
new dependences appeared in the process, and the stan- 
dard deviation of the CORAVEL-CfA minus standard 
radial-velocity differences decreased from 0.38 to 0.33 
km s for the stars in common. This quantity can be 



viewed as an approximate estimate of the systematic 
accuracy of the CORAVEL-CfA catalog. 

We used the same procedure to bring other cata- 
logs to the IAU.Keck standard system. As for the stars 
of the CORAVEL-CfA catalog, the most widespread 
systematic errors for these stars depend on radial ve- 
locity and B—V. Dependences on these two parameters 
are usually manifest in observations of radial velocities 
(see, e.g., [27]). 

After transforming all the radial velocities to the 
single system, we calculate for each star its weighted 
mean radial velocity. In this case, we set the weights in 
accordance with the standard deviation of the catalog 
minus standard radial-velocity differences for stars in 
common to the catalog and standard. We assign sin- 
gle observations the same weight as that for catalogs 
obtained from observations performed with the same 
instrument. 

The radial velocities of the CORAVEL-CfA catalog 
transformed to the standard IAU.Keck system have 
such high precisions compared to many other sources 
and the CORAVEL-CfA catalog has such high weights 
that the weighted mean velocities for these stars barely 
differ from the corresponding CORAVEL-CfA radial 
velocities with the systematic errors corrected. 

However, some of the CORAVEL-CfA data should 
be treated with caution. Figure 5 compares the radial 
velocities and U, V, W space- velocity components from 
the CORAVEL-CfA catalog and from our calculations 
based on other two sources of radial velocities for 3400 
stars in common (all radial velocities were transformed 
to a single system and the same HIPPARCOS/Tycho-2 
proper motions were used). All the quantities are ap- 
preciably correlated. The differences in the radial ve- 
locities are due mostly to errors in catalogs that are 
less accurate than the CORAVEL-CfA catalog. How- 
ever, a fairly large number of stars exhibit discrepances 
between the U, V, and W values from the CORAVEL- 
CfA catalog and our calculations. This discrepancy is 
due to the distances employed: we used distances calcu- 
lated from the HIPPARCOS trigonometric parallaxes, 
whereas photometric parallaxes were used to calculate 
the distances in the CORAVEL-CfA catalog and the 
corresponding factor for transforming the proper mo- 
tions into linear velocities for stars for which the rel- 
ative errors of the HIPPARCOS parallaxes is greater 
than 13%. The authors of the CORAVEL-CfA cata- 
log admit that these photometric parallaxes are not 
suitable for binary stars, giants, and several other cat- 
egories of stars. We found most of the stars with fairly 
discordant U, V, and W velocity components for the 
CORAVEL-CfA catalog and our own data, indeed, 
to be such "problematic"objects, whose photometric 
parallaxes in the CORAVEL-CfA catalog are system- 
atically greater than the corresponding trigonometric 
parallaxes (and whose heliocentric distances are sys- 
tematically smaller). As is evident from Fig. 6, the 
same effect arises for the same "problem "stars when 



we compare the distances based on the HIPPARCOS 
parallaxes and the CORAVEL-CfA distances, and in 
the appreciable difference between our absolute mag- 
nitudes and the absolute magnitudes adopted from the 
CORAVEL-CfA catalog (note that the discrepancy be- 
tween the interstellar extinctions inferred using differ- 
ent methods is insignificant, since we are dealing with 
fairly nearby stars whose extinction does not, on aver- 
age, exceed O.l" 1 ). 

Among the stars considered here are more than 
7000 known or suspected multiple stars, including more 
than 800 known and 4500 suspected spectral binaries. 
Therefore, one of the problems we had to face when 
calculating the weighted mean radial velocities were 
the errors in the identification of components of multi- 
ple systems made in many publications. In other words, 
authors often incorrectly identify the component whose 
radial velocity they have measured, or even fail to men- 
tion the component to which the measurement refers. 
These errors are largely due to the fact that, until 
recently, the nomenclature for stellar and planetary 
systems was not standardized; i.e., the same compo- 
nent had different designations in different catalogs. 
We made a number of corrections of this type. 

In all cases where there is no doubt that the compo- 
nents of a stellar system are gravitationally bound (the 
component separation does not exceed 20 000 AU), we 
calculated only the systemic radial velocities, and will 
consider only these velocities further. In the case of 
gravitationally wide stellar pairs (with component sep- 
arations exceeding 20 000 AU), the OS AC A catalog 
lists all components with measured radial velocities. 

Estimates of the errors in the radial velocities are 
available for most of the stars in the OSACA database, 
enabling us to estimate the random errors in the ob- 
served space velocities. For example, 26 492 stars sat- 
isfy the criteria 

e OT /7r < 0.3, (1) 
e v < 15 km s _1 , (2) 
V pec < 100 km s" 1 , (3) 

Here, e w is the fractional error in the parallax, ey = 

\J £y + z\ t + ey b is the error in the observed velocity 

vector V = (Vr, VI, Vb), where V r is the radial compo- 
nent of the velocity, V = 4.74^ cos br is the velocity 
component in the direction of Galactic longitude, and 
Vb = 4.74/i;,r is the velocity component in the direc- 
tion of Galactic latitude; V pec = \JU 2 + V 2 + W 2 is 
the residual (peculiar) velocity of the star relative to 
the local standard of rest. We estimated the errors in 
Vi and Vb taking into account the errors in the paral- 
lax but without allowance for correlations between the 
proper-motion components [iicosb and /ib, using the 
formula 



4.74 



Hv,,v b ) 



where are the tangential-velocity components. We 
set the velocity of the Sun relative to the local standard 
of rest equal to (%,%,%) = (10.0,5.3,7.2) km s -1 , 
in accordance with [28]. We found the mean random 
errors to be ey r = 2.7 km s _1 , ey, = 2.3 km s _1 , and 
ey b = 1-7 km s _1 . 



4. METHODS FOR THE KINEMATIC 
ANALYSIS 

In this paper, we analyze the kinematics of stars 
based on the OSACA catalog in an Ogorodnikov-Milne 
model. We analyze the residual velocities using the clas- 
sic method of determining the Schwarzschild ellipsoid 
and two-dimensional distributions of the U, V velocity 
components. 

4-1- The Ogorodnikov-Milne Model 

We used a linear Ogorodnikov-Milne model with 
the notation introduced by Clube [29, 30]. The ob- 
served velocity V(r) of a star with a heliocentric radius 
vector r can be described by the following vector equa- 
tion, up to first-order in the small quantity r/R <C 1: 



V(r) =u +Mr + V, 

where u is the velocity of the centroid of the stars 
relative to the Sun (—V®(uq,v@,w@)), V is the resid- 
ual velocity of the star (here we assume that the stel- 
lar residual velocities are randomly distributed), and 
M is the matrix (tensor) of the displacements, whose 
components are equal to the partial derivatives of 
11(1*1,1*2,1*3) with respect to r(ri,r2,T"3): 



{ dup 
[dr q 



(p, q= 1,2,3). 



The matrix M can be decomposed into the local defor- 
mation tensor M + and the local-rotation tensor M~ 
[31]: 



(see [32] for more details). The set of conditional equa- 
tions (5) .(7) contains 12 unknowns: three velocity com- 
ponents for Vq and nine components of M pg , which 
are estimated using the least-squares method. We then 
calculate the components of the deformation and rota- 
tion tensors using the already derived M pg values based 
on (4). Thus, the Oort constants can be calculated as 
follows: A = 0.5(Mi2 + M 21 ), B = 0.5(M 2 i - M 12 ), 
C = 0.5(Mn - M 22 ) h K = 0.5(Mn + M 22 ), and the 
vertex deviation l xy can be calculated from the for- 
mula taa2l xy = —C/A. The angular velocity of the 
Galactic rotation can be obtained from the formulas 
uj = B — A = —M\2. The form of these equations is 
such that negative oj corresponds to the direction of 
the Galactic rotation. 

4-2. Determination of the Schwarzschild Ellipsoid 

We calculated the elements of the residual- velocity 
ellipsoid (Schwarzschild ellipsoid) using the well- 
known statistical method described in detail in [31, 33, 
34]. This method consists in determining the symmetric 
momentum tensor. We have six equations for each star 
that can be used to determine the six unknown com- 
ponents of the velocity-dispersion tensor. Analysis of 
the eigenvalues of the velocity dispersion tensor yields 
the principal semiaxes 0X2,3 of the residual-velocity el- 
lipsoid (71,2,3 and also the directions ^1,2,3, &i,2,3 of the 
principal axes. When writing the formulas for the stel- 
lar residual velocities, we take into account the overall 
rotation of the Galaxy by adopting the Oort constants 
A = 13.7km s" 1 kpc" 1 and B = -12.9 km s" 1 kpc" 1 , 
in accordance with [35]. 

4-3. Two-dimensional Maps of the Space Velocities 

It is of interest to construct and analyze two- 
dimensional velocity maps, in particular, velocity maps 
in the U,V plane [36, 37], in order to analyze the 
fine structure of the distributions of the residual space 
velocities in order to identify clusters, streams, etc. 
We estimated the two-dimensional probability densi- 
ty f(U,V) from the calculated, discretely distributed 
U, V velocities using the method of adaptive smoothing 
[36] with a two-dimensional radially symmetric Gaus- 
sian kernel: 
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(p, 9=1,2,3). 



' q KJI v . 
The working equations have the form 

V r = f{u & ,VQ,W Q ,M n , ...,M 33 ), 

Vi = /(u ,« ,Mii,...,M 23 ) ) 
Vb = f(u Q ,VQ, wq, M11, ...,M 33 ), 



(4) 



(5) 
(6) 

(7) 



where r 2 = x 2 + y 2 . In this case, the necessary rela- 
tion J K(r)dr = 1 for probability density estimates is 
satisfied. The main idea behind this method is that, 
at each point of the map, smoothing is performed by 
a beam whose size is determined by a parameter cr, 
which varies depending on the data in the neighbor- 
hood of the point considered. Thus, the smoothing is 
performed using a fairly narrow beam in high-density 
regions, and the beam width increases with decreasing 
density. We use the following definition for adaptive 



smoothing at an arbitrary point £ = (U, V) [36]: 

2—1 V 7 

where & = (L/j, V^), Ai is a local dimensionless parame- 
ter of the beam at the point £, , h is a general smooth- 
ing parameter, and n is the number of data points 
& = (Ui,Vi). The parameter Aj at each point of the 
two-dimensional U, V plane is defined as follows: 

where g is the geometric mean /(£)■ It is evident that, 
to determine Ai, we must know the distribution /(£), 
which, in turn, can be determined if all Ai are known. 
Therefore, the problem of finding the unknown distri- 
bution must be solved iteratively. As a first approxi- 
mation, we use the distribution obtained by smoothing 
the initial U, V map with a fixed-sized beam. The op- 
timum value of the parameter h can be determined by 
minimizing the mean square deviation of the estimate 
/(£) from the true distribution /(£). The h values that 
we inferred for the various stellar samples considered 
here vary from 6 to 10 km s _1 . A typical uncertainty 
in the velocity is 2 km s -1 in this case. This influ- 
enced our choice of the discretization interval for the 
two-dimensional maps; the area of a square pixel was 
taken to be S = 2x2 km 2 s~ 2 . The map must be scaled 
using the factor nS. 



5. RESULTS OF THE KINEMATIC ANALYSIS 

5.1. Effect of Binary Stars 

The aim of this section is to establish how the bina- 
rity of stars affects the parameters of the Ogorodnikov- 
Milne model. To identify possible features that de- 
pend on stellar age, we selected from the OSACA cat- 
alog stars of luminosity classes V and IV, i.e., main- 
sequence stars satisfying the conditions 

e n /iT < 0.1, ey < 8 km s -1 , 

and also condition (3). Figure 7 shows the distribution 
of random errors of the V r ,Vi,Vb components of the 
observed space velocities as a function of distance for 
the resulting sample of 6276 stars. We find the mean 
random errors of the velocity components to be ey r = 
1.1 km s -1 , ey, = 1.3 km s _1 , ey b — 0.9 km s _1 . In 
this case, the mean error in the absolute space velocity 
is 1.9 km s -1 . The homogeneous nature of our data 
enables us to solve the initial equations (5)-(7) with 
unit weights. 

We first subdivided the entire sample into two parts 
depending on the color index, with the boundary at 
B — V = 0.5, then produced for each part subsam- 
ples consisting of either single or binary stars. Table 



3 gives our derived components of the solar veloci- 
ty, Vq(uq,Vq, Wq) relative to the corresponding stellar 
centroid, coordinates Lq and B & of the apex of the so- 
lar motion, elements M pq of the displacement tensor 
and roots \i.2,3 of the deformation tensor, principal 
axes cr 12 ,3 of the residual-velocity ellipsoid, and direc- 
tions ?i,2,3)&i,2,3 of the principal axes. As is evident 
from the upper part of the table, the components of 
the solar motion based on the samples of single and bi- 
nary stars do not differ significantly in any of the color 
intervals. 

The only elements of the displacement tensor to dif- 
fer significantly from zero are M\% for the three samples 
indicated in columns 2-4 and M\\ for the last sample of 
binary main-sequence red dwarfs (column 5). The Mn 
values derived for the samples of single and binary red 
dwarfs differ by AMn ■ r w 3 ± 1 km s _1 . 

Two roots of the deformation tensor, Ai,2, differ 
significantly from zero for the two single-star samples 
(columns 2 and 4), whereas the third root, A3, does not 
differ significantly from zero. Hence, we have a planar 
case where = K ± V 'A 2 + C 2 [38]. The closeness 
of the absolute values of the roots Ai and A2 is sugges- 
tive of rotation (without expansion/ contraction). 

The deformation tensor for the binary red dwarfs 
represents a linear case with a single significant root 
of Ai = 117 ± 29 km s _1 kpc -1 , and with the fairly 
high value K = 55 ± 19 km s -1 kpc -1 . In this case, 
the ellipse of the deformation tensor is highly elongated 
along the x axis, i.e., we have expansion along a single 
axis. 

As is evident from Table 3, the elements of the 
residual-velocity ellipsoids derived for the single and 
binary samples do not differ significantly in any of the 
color intervals. The elements of the residual- velocity el- 
lipsoid derived for the binary stars have larger random 
errors. 

The samples considered are most refined in terms 
of the errors of the heliocentric distances and velocity 
components of the stars. 

Our distributions of U, V (Fig. 8) show that the bi- 
nary stars are much more concentrated than the single 
stars (by about a factor of two) near the local peaks 
associated with the Hyades and Pleiades clusters. Note 
that the bottom contour and contour step in the figures 
shown are 10% of the maximum value. 

The specific features of the parameters of the 
Ogorodnikov-Milne model found by analyzing the sam- 
ples of single and binary red dwarfs, such as the statis- 
tical significance of Mn, are consistent with our results 
based on an analysis of a significantly smaller sample 
of nearby visual binaries [39]. Note that we have de- 
termined the kinematic parameters with much high- 
er precision in this current work. The distributions of 
the U, V velocities of both the single and multiple stars 
(Fig. 8), and of the red dwarfs (the oldest stars), in par- 
ticular, show that star formation in the Galactic disk 
is influenced by some persistent factor, which could be 



associated with the bar at the Galactic center [40,41] 
and/or the spiral structure of the Galaxy [42]. 

5.2. Negative K effect 

The aim of this section is to analyze stars exhibit- 
ing the strongest negative K effect. To this end, we 
selected only single stars that obey relations (l)-(3) 
and the criterion 7r < 10 milliarcseconds, which leaves 
us mostly with giants. All the selected stars are located 
at distances of 0.1-0.6 kpc. Our analysis of 9468 stars 
(r = 0.175 kpc) of mixed spectral composition yield- 
ed A = 12.3 ± 1.0km s -1 kpc -1 , B = -13.4 ± 1.0 
km s -1 kpc , C = —6.1 ± 1.0 km s -1 kpc -1 , and 
K = —4.1 ± 1.0 km s -1 kpc -1 . We found the vertex 
deviation to be I = 13° ± 2° and the angular velocity 
of the Galactic rotation to be u — — 25.7 ± 1.2 km s -1 
kpc -1 . The derived values of the Oort constants A, B, 
and 0j o are, on the whole, fairly consistent with the 
Galactic rotation parameters obtained in various other 
studies based on much more distant stars than those 
used here [35, 43, 44]. At the same time, we found K to 
differ significantly from zero, which is of considerable 
interest. A negative K effect in the Ogorodnikov-Milne 
model means that the stellar system considered is in a 
state of contraction. 

Several authors [35, 45] have found negative K ef- 
fects of -(1-6) km s -1 kpc -1 in the motion of OB stars 
(out to heliocentric distances of « 4 kpc) . The nature 
of this result is not entirely clear. Various hypotheses 
have been suggested: it could be an effect due to the 
specifics of the methods used to measure the stellar ra- 
dial velocities [46], or a result of the influence of the 
bar in the center of the Galaxy [46] , or the effect of the 
spiral structure [47,48]. Fernandez et al. [45] tried to 
eliminate the negative K effect by taking into account 
the spiral structure, but their attempt failed, and the 
negative K effect remained. 

Note that the youngest and nearest OB stars exhib- 
it a characteristic positive K effect [32, 49], and there- 
fore we do not consider these stars here. 

In this work, we succeeded in forming a sample of 
stars with a maximum negative K effect, which consists 
of 1269 stars with color indices B - V < 0.2. We did 
not use stars of luminosity classes V and IV or O and B 
stars, so that the sample in question contains only stars 
of spectral types A0.A5 and luminosity classes I, II, III. 
However, the bulk of the sample is made up of stars 
with no luminosity classes given in the OSACA catalog. 
We then analyzed the residual velocities of the stars 
corrected for the overall Galactic rotation (see Section 
3.2). Our analysis yielded for the peculiar velocity of 
the Sunu© = 10.2±0.4 km s -1 , v Q = 10.9±0.4 km s -1 , 
w Q = 6.6 ±0.4 km s -1 , V Q = 16.4±0.4 km s -1 , L Q = 
47 ± 2°, B & = 24 ± 2°, and for the components M+ 
and M~ of the tensor M (in units of km s -1 kpc -1 ) 



characterizing the residual kinematic effects 
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The roots of the deformation tensor (Ai,A2,As) are 
(3.5,-29.3,-3.5) km s -1 kpc -1 , with standard errors of 
-3 km s -1 kpc -1 . The directions of the principal axes 
of the deformation tensor are 

Li = 282 ±2°, Bi=-12±l°, 
L 2 = 12±4°, B 2 = -l±26°, 
L 3 = 285 ±4°, B 3 = 79±6°. 



The vertex deviation in the xy plane is l xy = —33 ±4°; 
we also derived the parameters of the K effect: K = 
0.5(Mn + M 22 ) = -13.1 ±2.0 km s -1 kpc -1 , and C = 
0.5(Mn -M 22 ) = -14.9 ±2.0 km s -1 kpc -1 . 

The rotation tensor, M~ , contains no elements that 
differ significantly from zero, leading us to conclude 
that the radial component is the dominant effect in the 
residual motions of the stars considered. In this case, 
we set one of the roots (A3) equal to zero and consider 
only the xy plane. Since ((dVR,/dR)R = Ai, we have 
(Vr/R)r = A2, where R is the unknown distance from 
the kinematic center. However, the relations 

if + C = 4±3kms -1 kpc -1 , 
K - C = -29 ± 3 km s -1 kpc -1 , 

yield an estimate of the characteristic residual velocity 
of the stars, (K — C) ■ r = —5.4 ± 0.6 km s , and 
also K = —12.9 km s -1 kpc -1 and C = 16.4 km s -1 
kpc -1 (the sign has changed), which characterizes the 
radial component of the residual stellar velocities in the 
case considered. The resulting estimate, —5.4 ±0.6 km 
s -1 , is consistent with the amplitude of the periodic 
radial oscillations of the stellar residual velocities in 
the interarm space, Jr = 7 ± 1 1 km s -1 , found by 
Mel'nik et al. [50] based on an analysis of the motions 
of OB associations. 

The direction of the third axis, B3 — 79 ± 6°, indi- 
cates compression in the Galactic plane. In fact, both 
roots of the deformation tensor (Ai and A 2 ) are close 
to zero. This means that the compression occurs along 
a single axis in the direction 12 — 192°, with its error 
being 2°-4°. 

The principal axes of the residual-velocity tensor 
(0-1,0-2,0-3) are (16.09 ± 0.56, 10.04± 0.53, 6.89 ± 0.47)) 
km s -1 , and the directions of the principal axes have 



the coordinates 

Zi = 27±4°, 6i = -l±0°, 
/ 2 = H7±6°, b 2 = -l±2°, 
Z 3 = 80±6°, 6 3 = 88 ±3°. 

We can see that the angles h, l xy and L2 differ from 
each other, which further emphasizes the fact that the 
effect of the radial component is significant in the kine- 
matics of these stars. 

Figure 8c shows the distribution of the residual U, V 
velocities for the stars considered. The stars are con- 
centrated in only two peaks. the so-called Sirius and 
Pleiades peaks. It is striking that the velocity of the 
stars displaying a negative K effect relative to the lo- 
cal standard of rest (vq) is fairly small. 

Rybka's [51] analysis of the proper motions of about 
40 000 giants of luminosity class III yielded K = 
—7.5 ± 1.8 km s _1 kpc -1 for G5-K0 stars located be- 
yond 250 pc (with their mean distance being about 
364 pc). Our result agrees well with that of Rybka [51], 
leading us to conclude that the use of the high-precision 
space velocities of the OSACA stars has enabled us to 
detect this negative K effect at much closer distances. 
Our approach has also shown that we actually have 
compression along a single axis, enabling us to come 
closer to understanding of the nature of this effect. 



6. CONCLUSIONS 

We have created and have been continuously up- 
dating the OSACA (Orion Spiral Arm CAtalogue) 
database for stars with known coordinates, parallax- 
es, proper motions, and radial velocities located in the 
Orion arm of our Galaxy. We transformed radial ve- 
locities adopted from over 1400 sources, including the 
Geneva-Copenhagen survey of the solar neighborhood 
(CORAVEL-CfA), into a single system of radial veloc- 
ities based on 854 standard stars from our list. This 
enabled us to calculate the weighted mean radial ve- 
locities of more than 25 000 stars of the HIPPARCOS 
catalog with a median error of 1 km s -1 . We use the 
elimination of systematic errors of the CORAVEL-CfA 
catalog as an example to demonstrate the efficiency of 
the proposed standard of stellar radial velocities. 

We have corrected numerous errors in the identi- 
fication of components of multiple systems made in 
many publications. We have used the OSACA data to 
analyze the kinematics of 6276 nearby main-sequence 
stars with space- velocity errors ±2 km s" 1 . The use of 
an Ogorodnikov-Milne model and analysis of the dis- 
tribution of the U, V velocities the of stars revealed 
significant differences in the kinematics of single and 
multiple stars. In particular, multiple stars are much 
more concentrated than single stars in the vicinity of 
the peak associated with the Hyades-Pleiades stream 
in the U, V velocity plane. 



Based on a sample of 9468 relatively distant 
((r w 0.2 kpc) stars of various spectral types, we 
estimated the angular velocity of the Galactic rotation 
to be LUo = — 25.7 ± 1.2 km s _1 kpc -1 (we found the 
vertex deviation to be I — 13 ± 2°). The radial compo- 
nent of this velocity indicates a K effect, which shows 
up most conspicuously in the motions of A0-A5 gi- 
ants. We found for the parameters of the Ogorodnikov- 
Milne model K = -13.1 ± 2.0 km s" 1 kpc" 1 and 
C = 14.9 ± 2.0 km s" 1 kpc" 1 , which we used to es- 
timate the characteristic residual velocity of the stellar 
sample considered, —5.4 ±0.6 km s . The maximum 
compression is in the direction 12 — 192°, with an error 
of 2° -4°. 
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Table 1. Distribution of intrinsically bright stars (in percent) in the Galactic octants and heliocentric distance 
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13.3 
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813 





Fig. 2. Distributions of OSACA stars in absolute magnitude and heliocentric velocity, for single (solid) and 

multiple (dashed) stars. 



Table 2. Format of the current release of the OSACA catalog. 
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Table 3. Kinematic parameters of nearby main-sequence stars. 
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Fig. 4. CORAVEL-CfA minus standard radial-velocity differences for 484 stars in common as a function 
various parameters before (left-hand plots) and after (righthand plots) allowance for systematic effects. 
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Fig. 5. Comparison of the (a) radial velocities, (b) U velocity components, (c) V velocity components, and (d) 
W velocity components from the CORAVEL-CfA catalog with the corresponding velocities calculated from 
other radial-velocity sources (a total of 3400 stars in common) for. 
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Fig. 6. Left: comparison of absolute magnitudes from the CORAVEL-CfA catalog with the absolute 
magnitudes calculated from the HIPPARCOS parallaxes with allowance for interstellar extinction (3400 stars 
in common) . The difference is due to the differences in the allowance for extinction and the use of inaccurate 
photometric parallaxes in the CORAVEL-CfA catalog. Right: heliocentric distances of stars of the 
CORAVEL-CfA catalog compared to their heliocentric distances from the HIPPARCOS catalog. 




Fig. 7. Random errors of the velocities of nearby main-sequence stars as a function of heliocentric distance. 
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Fig. 8. Distributions of U and F (relative to the local standard of rest) for (a) single stars, (b) nearby main 
sequence binaries with (B — V > 0.5), and (c) stars with a negative K effect. 



